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Abstract 

We calculate the contributions of the Kaluza-Klein (KK) modes to the 7-penguins, gluon- 
penguins, 7-magnetic penguins and chromomagnetic penguins in the Appelquist, Cheng and 
Dobrescu (ACD) model with one universal extra dimension. Together with our previous calcu- 
lation of penguin diagrams [1] this allows to study the impact of the KK modes on the decays 
B As7, B ^ Xg gluon, B XgH^jj.^ and Ki 7r^e~^e~ and on the CP-violating ratio e' /e. 
For the compactification scale 1 /R = 300 GeV the perturbative part of the branching ratio for 
B XsfJ,'^ iJ-~ is enhanced by 12% while the zero in the ^fb asymmetry is shifted from sq = 0.162 
to So = 0.142. The sizable suppressions of Br{B Xsj) (20%) and Br{B Xs gluon) (40%) 
could have interesting phenomenological implications on the lower bound on l/i? provided the 
experimental and theoretical uncertainties will be decreased. Similar comments apply to e' /e 
that is suppressed relative to the Standard Model expectations with the size of the suppression 
depending sensitively on the hadronic matrix elements. The impact on Kl — > vr'^e+e" is below 
10%. We point out a correlation between the zero sq in the Ayb asymmetry and Br{B Xg^) 
that should be valid in most models with minimal flavour violation. 



1 Introduction 



In a recent paper [1] we have calculated the contributions of the Kaluza-Klein (KK) modes to 
the Kl — Kg mass difference AM/^-, the parameter ek-, the — ^ mixing mass differences 
AMd^s and rare decays tt^vu, Ki^ ir^uu, Kj^ fi^fi~, B Xs^cfi^u and Bg^d 

in the Appelquist, Cheng and Dobrcscu (ACD) model [2] with one universal extra dimension. 
In this model the only additional free parameter relative to the Standard Model (SM) is the 
compactification scale 1/R. Thus all the masses of the KK particles and their interactions among 
themselves and with the SM particles are described in terms of 1/R and the parameters of the 
SM. This economy in new parameters should be contrasted with supersymmetric theories and 
models with an extended Higgs sector. 

A very important property of the ACD model is the conservation of KK parity that implies 
the absence of tree level KK contributions to low energy processes taking place at scales fi <^ 1/R. 
In this context the flavour changing neutral current (FCNC) processes considered in [1] are of 
particular interest as they appear in the SM first at one-loop and are strongly suppressed. 
Consequently the one-loop contributions from the KK modes to them could in principle be 
important. Indeed our analysis in [1] has shown that for < 400 GeV the KK effects in some 
of the processes in question could be significant. This is in particular the case of the decays 
governed by the CKM element \Vts\, like B — Xgi'u and Bg — > which the enhancement 

through the KK modes in penguin diagrams is not softened by the suppression of the relevant 
CKM parameters in contrast to the processes governed by \Vtd\- 

In the present paper we extend the analysis of [1] by calculating the decays that receive 
contributions from the 7 penguins, gluon penguins, 7-magnetic penguins and chromomagnetic 
penguins. In particular we analyze in detail the impact of the KK contributions on the decays 
B Xgj, B ^ Xg gluon, B XsH'^n~ and Kl 'K^e^e~ and on the CP-violating ratio 
e' /e. Among these processes only the decay B Xg^y has been so far considered in the ACD 
model in the literature [3] . As we will see our result is consistent with the one obtained by these 
authors but differs in details. 

Of particular interest are the decays B Xgj and B XgiJ.'^/j." as the hadronic uncer- 
tainties in these decays are moderate and the experimental branching ratios are found to be in 
the ballpark of the SM expectations. Thus in principle some constraints on the compactifica- 
tion scale 1/R could be obtained from these decays. On the other hand the upper bound on 
Kl — -K^e^e' is still two orders of magnitude above the predictions of the SM and the decay 
B ^ Xg gluon and the ratio e'/e are subject to very large hadronic uncertainties. Still, as we 
will see below, our results for B ^ Xg gluon and e' /e could turn out to be important in the 
future. 
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The effects of the KK modes on other processes of interest have been investigated in a 
number of papers. In [2, 4] their impact on the precision electroweak observables assuming a 
hght Higgs {niH < 250 GeV) and a heavy Higgs led to the lower bound l/R > 300 GeV and 
1/R> 250 GeV, respectively. Vacuum stability in a simplified ACD model has been examined 
in [5]. 

Subsequent analyses of the anomalous magnetic moment [6], and of the Z ^ bb vertex [7] 
have shown the consistency of the ACD model with the data for l/R > 300 GeV. The latter 
calculation has been confirmed by us [1] . 

The scale of l/R as low as 300 GeV would lead to an exciting phenomenology in the next 
generation of colliders [8, 9, 10, 11, 12]. Moreover the cosmic relic density of the lightest KK 
particle as a dark matter candidate turned out to be of the right order of magnitude [13]. The 
related experimental signatures have been investigated in [14]. In the present paper, as in [1], 
in order to be more general we will include the results for 1/R = 200 GcV that is only slightly 
below the lowest value of l/R = 250 GcV allowed by electroweak precision data. 

As our analysis of [1] shows, the ACD model with one extra dimension has a number of 
interesting properties from the point of view of FCNC processes discussed here. These are: 

• GIM mechanism [15] that improves significantly the convergence of the sum over the KK 
modes corresponding to the top quark, removing simultaneously to an excellent accuracy 
the contributions of the KK modes corresponding to lighter quarks and leptons. This 
feature removes the sensitivity of the calculated branching ratios to the scale Mg ^ l/R 
at which the higher dimensional theory becomes non-perturbative and at which the towers 
of the KK particles must be cut off in an appropriate way. This should be contrasted with 
models with fermions localized on the brane, in which the KK parity is not conserved and 
the sum over the KK modes diverges. In these models the results are sensitive to Mg and 
the KK effects in AMg^d are significantly larger [16] than found in [1]. We expect similar 
behaviour in the processes considered here. 

• The low energy effective Hamiltonians are governed by local operators already present in 
the SM. As flavour violation and CP violation in this model is entirely governed by the 
CKM matrix, the ACD model belongs to the class of the so-called models with minimal 
flavour violation (MFV) as defined in [17]. This has automatically the following important 
consequence for the FCNC processes considered in [1] and in the present paper: the impact 
of the KK modes on the processes in question amounts only to the modification of the 
Inami-Lim one- loop functions [18]. 

• Thus in the case of the processes considered in [1] these modifications have to be made 
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in the function S [19] in the case of AM^^g and of the parameter ek and in the functions 
X and Y [20] in the case of the rare decays analyzed there. In order to study the decays 
B — Xgj, B ^ Xg gluon, B — XgH'^n" and Kl — Tr^e'^e" and the CP-violating ratio 
s'/e the KK contributions to new short distance functions have to be computed. These 
are the functions D (the 7 penguins), E (gluon penguins), D' (7-magnetic penguins) and 
E' (chromomagnetic penguins). In the ACD model these functions depend only on mt 
and the single new parameter, the compactification radius R. 

Our paper is organized as follows. In section 2, we summarize those ingredients of the ACD 
model that are relevant for our analysis. In particular, we give in appendix B the set of the 
relevant Feynman rules involving photons and gluons that have not been given so far in the 
literature. Further details can be found in [1, 2]. In section 3, we calculate the KK contributions 
to the penguin functions D, E, D' and E'. In section 4, we analyze the decay B Xgj and 
compare our result with the one of [3]. In section 5, we discuss the decay B ^ Xg gluon. The 
corresponding analyses of 5 — XgjJL^iJr , Kl — 7r°e"^e~ and e'/e are presented in sections 6, 
7 and 8, respectively. In section 9, we summarize our results and give a brief outlook. As a 
byproduct of our work we point out a correlation between the zero sq in the ^fb asymmetry 
in B ^ XsiS^jj," and Br{B Xgj) that should be valid in most models with minimal flavour 
violation. We also generalize the background field method to five dimensions. 

2 The Five Dimensional ACD Model 

The five dimensional ACD model [2] has been described already in detail in our previous paper 
[1] where we have given all Feynman rules necessary for the calculations considered there. There 
are two new features in this work. The first one is the addition of QCD in 5 dimensions which is 
straightforward to formulate and has already been discussed in some detail in [21]. The other one 
is the background field method described in appendix A which is used in the off-shell matching 
calculation. The additional Feynman rules that we need in the present paper, namely those 
involving photons and gluons, are collected in appendix B. 

Here we only recall the most important features of the ACD model. The fifth dimension is 
compactified on the orbifold S^/Z2 to reproduce chiral fermions in 4 dimensions. In addition 
to the ordinary particles of the SM, denoted as zero modes (n = 0), there are infinite towers of 
KK modes (n > 1). There is one such tower for each SM boson and two for each SM fermion, 
while there also exist physical scalars 

«(n) 

and a^^j with (n > 1) that do not have any zero mode 

partners. The masses of the KK particles are universally given by 
2 

mfn)=ml + ^. (2.1) 
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Here mo is the mass of the zero mode, as Mw, Mz, mt- For a°^^ and this is Mz and M\y, 
respectively. The compactification radius R is the only additional parameter relative to the SM. 
To a quark SU{2)l doublet, there correspond the KK modes Q, whereas the partners of the 
right-handed singlets are U and V. 

As the KK modes contribute to the process considered here first at the one loop level, the 
natural variables that enter the Inami-Lim functions are 

i{n) 1^2.2) 



W(n) 

with denoting the masses of the fermionic KK modes and M^(„) the masses of the W 

boson KK modes. However, in phenomenological applications it is more useful to work with the 
variables xt and Xn defined through 

_ m'f _ rrin _ n 

than with 



3 Penguin Diagrams in the ACD Model 

3.1 Preliminaries 

The rare decays in the ACD model considered here are governed as in the SM by various penguin 
diagrams. The SM contributions to AF = 1 box diagrams are subleading but non-negligible. 
On the other hand the KK contributions to the latter diagrams are tiny [1] and can be neglected. 

The penguin vertices including electroweak counter terms can be conveniently expressed in 
terms of the functions C (Z^ penguins), D (7 penguins), E (gluon penguins), D' (7-magnetic 
penguins) and E' (chromomagnetic penguins). In the 't Hooft-Feynman gauge for the and 



propagators they are given as follows: 

s^d = - i\^^D{xt, l/i?)s-((z'7M - % - (3.2) 

sG'^d = - iX^^^Eixt, l/R)-Sa{<f^>^ - % m - l^)T^pd0 (3.3) 

s^h = zAt^^L»'(xt,l/i?)sKA9^K(l+75)]]6 (3.4) 

sG'% = i\t%-^E\xul/R)U^(^^xq\m^{l+^^W^php, (3.5) 
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where is the Fermi constant, Oy, is the weak mixing angle and 



(3.6) 



In these vertices is the outgoing gluon or photon momentum and T" are colour matrices. The 
last two vertices involve on-shell photon and gluon, respectively. We have set = in these 
vertices. 

Each function in (3.1)-(3.5) is given by 



F{xul/R) = Fo{xt) + Y,Fn{xuXn), F = C,D,E,D',E', 

n=l 

with Xn defined in (2.3). The functions F„(x(,x„) in (3.7) are defined through 

Fnixt,Xn) = G(xt(n)) - G(x„(„)), 



(3.7) 



(3.8) 



with a:j(„) defined in (2.2) and the functions G(x((„)) and G(x„(„)) representing the contributions 
of the Qt{n), Ut{n) and Qu{n)^ ^u{n) modes, respectively. 

The functions Fo{xt) result from the penguin diagrams in the SM and the sum represents the 
KK contributions to the relevant penguin diagrams. Fo{xt) were calculated by various authors, 
in particular by Inami and Lim [18]. They are given explicitly as follows: 



Coixt) 



Xt 



xt-6 3xt + 2 



_ , , 4 , -19a;? + 25a;? a;?(5a;? - 2xt - 6) , 

, 2, x?(15- 16xt + 4a;?) , xt(18 - llx* - 

Eoixt) = -- Inx. + 6(l-x.)^ + 12(1 -x,)3 



D'oixt) 



(8xf + 5xg - 7xt) xg(2-3xt) 
12(1 -xt)3 + 2(l-xt)4 

xt{xf — 5xt — 2) 3 xf 



+ 



Inxt . 



4(1 - xt)3 2 (1 - xt)4 
The Z° penguin functions Cn{xt,Xn) have been calculated in [1] with the result 



Cfi {xt , Xji ) — 



Xt 



8{xt - 1)2 



Xf — 8xt + 7 + (3 + 3x( + 7xji — xtx^ )ln ' " 

i+Xn 



(3.9) 

(3.10) 
(3.11) 
(3.12) 
(3.13) 

(3.14) 



In the present paper we will calculate the remaining functions Dn{xt, x^), Enixt, Xn), D'^{xt, Xr. 
and E'^{xt,Xn). 
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3.2 General Structure of the Calculation 



The function Fn{xt, x„) with F — D,E, D' , E' can be found by calculating the vertex diagrams 
in fig. 1. Contrary to the penguins where one has to add an electroweak counter term as 
discussed in [22], this is not necessary for the 7 and gluon penguins. Here the counter terms 
are only formally used to render zero the coefficients of the dimension 4 operators sj^PlQ and 
sT^-^PlQ with q = d,b. This is a consequence of gauge invariance when the quark fields are 
set on-shell with the renormalization condition given in [23]. 






\ 



(4) 



(5) 



(6) 



Figure 1: Penguin diagrams contributing to Fn{xt,Xn)- 

In contrast to the calculation of the Z^-vertex the external momenta in fig. 1 and the masses 
of external quarks cannot be neglected. 

The calculation of the functions F was done in two different ways. The common starting 
point is the effective Hamiltonian at the matching scale nw 

' 36 

Ci{lJ,w)Qi + C7^il^w)Q7'y + C8Gil^w)Q8G , (3.15) 



nes{b^sj,G) = --^v;,Vtb 

with the operators [24] 



.1=30 



Qvj = ^mb{sa^.PRb)F^'-', 



(3.16) 
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Q8G = 






(3.17) 


Q30 = 






(3.18) 


Q31 = 






(3.19) 


Q32 = 








(3.20) 


Q33 = 






(3.21) 


Q34. = 






bG'''^"', (3.22) 


Q35 = 






bPf", (3.23) 


Q36 = 


^{sPRJi,b)d,F'''' - Qg. 




(3.24) 



The covariant derivative in the effective Lagrangian is defined as^ 
^/.V' = {d^^ + ieQ^A^ + igsT'^GD^. (3.25) 

In the operators Q34 and Q35 it acts only on the spinors, but not on the field strength tensors 
F^''' and G''^^"' . 

In this off-shell operator basis we have omitted the four-quark operators as they are not 
relevant for the calculation of the functions F. They can be found in [24], in particular the 
operators Q4 and Qg that appear in the definition of Q31 and Qz&- The effective vertices (3.2)- 
(3.5) correspond to the operators Qae, Q31, Q77 and QsGj respectively. 

In the first method the equations of motion of the quark fields and the relation = are ap- 
phed both on the full and on the effective side of the matching amplitude. After renormahzation 
the amplitude on the effective side is 

AT' = [^'(2^1 + T2) - DT,] b (3.26) 

with Ti = rrib^^f^PR and T2 = -2mbqij,PR. 

The full side was evaluated with the Feynman diagrams (l)-(6) in fig. 1 in 't Hooft-Feynman 
gauge. The vertices v4(o)W^(^)fl(^) ^^id ^(o)G^)a^^ are zero, hence there are no diagrams with 
these vertices to be considered. The result can be written as 



-4fu7r = [^1^1 + ^2^2] b, (3.27) 

yielding 



D = ^Hi-H2, D' = ^Hi. (3.28) 



^On the effective side we use a different sign convention for the couphngs e and gs compared to the Feynman 
rules in this paper and in [1]. 
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For the functions E and E' the procedure is analogous. As the gluon couples only to the 
quarks the relevant diagrams in this case are (1) and (2). 

The second method uses the off-shell amplitude for the matching, i.e. no equations of motion 
are applied. In order to maintain gauge invariance the method of background fields [25] is used. 
Here the matching is done with the full set of operators (3.16)-(3.24). 

For the full side the diagrams (l)-(4) in fig. 1 were evaluated in 't Hooft-Feynman gauge 
with the Feynman rules for the background fields 7 and G given in appendix B. As a virtue of 
the background field (BF) method, the diagrams (5) and (6) do not show up because the vertex 
^W^H^^fn) vanishes. 

The effective side can be directly matched to the full side. Comparing (3.2)- (3.5) with (3.15) 
we can read off 

-(0)^ ^ T^'-on(o)i 



D = -Ct^im), D' = ^2C'^'^{^iw). (3.29) 



E = -Cf^ ifiw), E' = ~2d'!l ii^w), (3.30) 



where "(0)" indicates the leading coefficients without QCD corrections. We found the same 
results for the functions F with both methods. 

3.3 The functions D, E D' E' and Z 

Adding up the contributions from the diagrams of fig. 1 we find 

_ (35 + 8a;t-19xf + 6a;2 (l0-9xf + 3x?) +3x„(53 - 58xt + 21x^)) 
" 108(a;, - if 

+ ^-[4-2xn + 4xl + xl) In-^"- 



+ Xji 

-4 (12 - 38xt + 54a;f - 27^ + 3x^ + 4 (3 + xt) + 3x^ (4 - + x?) 



En{xttXn) — 



18(xt - 1) 

A2xt-33x'} + 9x^]]\n — 

-\- Xfi 

xt (35 + 8xt - 19xf + 6x1 (lO - 9^t + 3a;?) + Sxn (53 - 58xt + 21x|)) 



+ Xn(-6 + 42xt - 33x| + 9x1) ) In ^!i±^ , (3.3I) 

^ ^ 1 -I- x„. 



36(xt - 1)' 
^ (1 + xn) (-2 + 3x„ + xl) In 



+ Xr. 



(1 + x„) (-6 + 19xt - 9x1 + x2 (3 + Xi ) + x„ (9 - 4xt + 3x|) ) x„ + xt 

H -. ^ ^In- , (3.32) 

6(xt - 1)^ 1 + x„ ' ^ ^ 
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xt (-37 + 4Axt + 17xf + 6x1 (lO - 9^;* + 3x|) - Sxn (21 - 54xt + 17x|)) 



, Xfi (2 IXji + ^Xjj) Xn, 

H :^ in ■ 



6 l + Xn 

(-2 + + 3xt) (xj + ?,x1 + Xn (3 + Xt) - X„ (1 + (-10 + Xt) Xt)) Xn + Xt 

m ■ 



6(-l + xt) 



l + xt 



, (3.33) 



E'n^XuXn) = 



Xt (-17 - 8xt + xl - 3xn (21 - 6xt + xj) - 6x^ (lO - 9xt + 3xf)) 

12(xt - 1)' 

--x„ (1 + x„) (-1 + 3xn) In " 
2 1 + a:„ 



(1 + Xn) (xt + 3xf +X^ (3 + Xt) - Xn (1 + (-10 + Xf)Xf)) X„ + Xt 

H ^In- . (3.34) 

2(xt-l)^ ^ + Xn ^ ' 

In fig. 2 we plot the functions F{xt, 1/R) versus 1/R. Tlie impaet of the KK modes on the 
function D is negligible. The function E is moderately enhanced but this enhancement plays 
only a marginal role in our phenomcnological applications. The most interesting are very strong 
suppressions of D' and E' , that for 1/R = 300 GeV amount to 36% and 66% relative to the SM 
values, respectively. 

While the AF = 2 box function S and the penguin functions E, D' and E' arc gauge 
independent, this is not the case for C, D and the AF = 1 box diagram functions B^'^ and B^^^ 
considered in [1]. In phenomenological applications it is more convenient to work with gauge 
independent functions [20] 



X(xt, l/R) = C{xt, l/R) + 5^^(xt, l/R) = Xo(xt) + AX 

y(xt, l/R) = C(xt, l/R) + B^f^ixt, l/R) = Yo{xt) + AY , 
Z{xt, l/R) = C{xt, l/R) + ^D{xt, l/R) = Zo{xt) + AZ 



with (mt = 167 GeV) 



Xoixt) 



Xt 



Xt + 2 3xt — 6 



+ 



Xt - 1 (xt - 1) 



In Xf 



Yoixt) = ? 



Xt -4 



+ 



3xt 



Xt - 1 (xt-l 



Inxt 



1.526 , 



0.980 



(3.35) 

(3.36) 
(3.37) 

(3.38) 
(3.39) 



^ , , 1 , 18xf - 163xf + 259x? - 108xt 
Zo{xt) = --lnxt + ^ TT^^ ^ + 



144(xt - 1)3 
32x^ - 38xf - 15x2 + 18xt 



72(xt - 1) 



Inxf = 0.679 



(3.40) 
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(c) (d) 

Figure 2: The functions F{xt,l/R) and Fo(xt) with F = D,E, D' , E' plotted versus l/R. The constant 
dashed hues are the SM values. 



summarizing the SM contributions and AX, AF and AZ representing the corrections due to 
KK modes. Our analysis of [1] and the negligible KK contributions to the function D found 
here imply to an excellent approximation 

oo 

AX = Ay = AZ = ^C7„(xt,x„). (3.41) 

n=l 

In table 1 we give the values of the gauge independent functions X,Y calculated in [1] and 
of Z, D' , E' calculated here for different \/R and rrit = 167 GeV. Our results for the function 
S have been confirmed in [26] . For completeness we give also the values of C and D in 't Hooft- 
Feynman gauge. We observe that for l/R = 300 GeV, the functions X, Y, Z are enhanced 
by 10%, 15% and 23% relative to the SM values, respectively. While these enhancements are 
smaller than the corresponding suppressions of D' and E' , the effect of the latter suppressions 
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will be softened in the relevant branching ratios through sizable additive QCD corrections. 



l/R [GeV] 


S 


X 


Y 


Z 


E 


D' 


E' 


C 


D 


200 


2.813 


1.826 


1.281 


0.990 


0.342 


0.113 


-0.053 


1.099 


-0.479 


250 


2.664 


1.731 


1.185 


0.893 


0.327 


0.191 


0.019 


1.003 


-0.470 


300 


2.582 


1.674 


1.128 


0.835 


0.315 


0.242 


0.065 


0.946 


-0.468 


400 


2.500 


1.613 


1.067 


0.771 


0.298 


0.297 


0.115 


0.885 


-0.469 


SM 


2.398 


1.526 


0.980 


0.679 


0.268 


0.380 


0.191 


0.798 


-0.476 



Table 1: Values for the functions S, X, Y, Z, E, D', E', C and D. 



4 B^Xsj 
4.1 PreliminEiries 

The inclusive B — > Xgj decay has been the subject of very intensive theoretical and experimental 
studies during the last 15 years. On the experimental side the world average resulting from the 
data by CLEO, ALEPH, BaBar and Belle reads [27] 

Br{B ^ Xsj)e,>im GeV = (3.281°:^^) • 10"^ . (4.1) 

It agrees well with the SM result [28, 29] 

Br{B ^ X,7)t;'>i.6 GeV = (3.57 ± 0.30) • 10"^ . (4.2) 

The most recent reviews summarizing the theoretical status can be found in [30, 31, 32]. The 
purpose of this section is the calculation of the impact of the KK modes on the branching ratio 
in question. 

As the NLO QCD corrections to Br{B Xsj) within the SM arc very important, it is 
mandatory to include them also in the ACD model. Fortunately, the QCD corrections for scales 
H < l/R are in the ACD model precisely the same as in the SM and the only elements of 
the NLO QCD corrections to Br{B Xgj) within the ACD model that arc missing, are the 
0{as) corrections to the functions D' and E'. The experience with such corrections within the 
SM and MSSM at low tan /? makes us to expect that they are significantly smaller than the 
renormalization group effects for fi < M^j\i and the 0{as) corrections at scales 0{mb). Therefore 
we arc confident that it is a good approximation to neglect these corrections in the case of the 
KK contributions. 
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Now, the NLO formulae for the relevant Wilson coefficients and the Br{B Xg^) are 
very complicated and not transparent. Therefore we will first discuss the impact of the KK 
contributions at the LO level. Subsequently we will present the numerical results of an NLO 
analysis that uses the formulae of [28] modified by the KK contributions to the functions D' 
and E' calculated in the previous section. 



4.2 Effective Hamiltonian 

The effective Hamiltonian for B — Xg'y at scales /x^ = 0{mb) is given by 



X] Ci{lJ.b)Qi + C7j{lJ.b)Q7j + CsG{lJ'b)QsG 



.1=1 



(4.3) 



where in view of | V*gVub |V^gVtb \< 0.02 we have neglected the term proportional to V*gVub- 
Here Q\....Qq are the usual four-quark operators whose explicit form can be found in [33]. The 
remaining two operators, characteristic for this decay, are the magnetic-penguins 

Q7i = ^mbsao^^'^il + i^)haF^,, QsG = |^m(,Saa^"(l + l^W^ph^G''^, (4.4) 

originating in the diagrams of fig. 1 with on-shell photon and gluon, respectively. 
The coefficients Ci(/Lt(,) in (4.3) can be calculated by using 

C{iib) = Ub{lJ-b,IJ^w)C{iiw) (4.5) 

Here U5{fib, fJ-w) with jjw = 0{Mw) is the 8x8 evolution matrix that is known including NLO 
corrections [34]. As stated before, we will first estimate the effect of the KK modes on C7^(/ife) 
and Csail^b) in the leading logarithmic (LO) approximation. In this case generalizing the SM 
formulae in [35] to the ACD model we obtain 

4f'^(M6) = r;i4°)(;x^) + ^(r?i-,i)c(g(/.^) + cf)(^H^)j:/.,ry-, (4.6) 

i=l 

8 

C^f"{f^b) = r?iC7^g(/XH^) + cf (/XM^)^M"% (4.7) 



i=l 



with 

asil^b) 



V = ^#4, (4.8) 



and 



(/XH^) = 1, c!f^M = -\D'{xt, l/R), C^%w) = -\e'{xu 1/R) (4.9) 
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and all remaining coefficients equal zero at /x = /j^w- The superscript "0" indicates the LO 
approximation. The "effective" coefficients in (4.6) and (4.7) are introduced in place of C7-y(/Xb) 
and CscilJ-b) in order to keep the LO Wilson coefficients renormalization scheme independent 
as discussed in [35]. The functions D'{xt, l/R) and E'{xt, 1/-R) have been calculated in section 
3. Finally the values of Oj, hi and hi are given in table 2. 



i 


1 


2 


3 


4 


5 


6 


7 


8 


tti 
hi 
hi 


14 

23 

2.2996 
0.8623 


16 
23 

-1.0880 



6 
23 
3 
7 




12 
23 
1 

14 




0.4086 
-0.6494 
-0.9135 


-0.4230 
-0.0380 
0.0873 


-0.8994 
-0.0185 
-0.0571 


0.1456 
-0.0057 
0.0209 



Table 2: Magic Numbers. 



Using the leading ^Ufe-dependence of ag in an effective five quark theory 
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and setting a^^\Mz) = 0.118 we find the results in table 3. 



(4.10) 





fib = 2.5 GeV 


/Xfe = 5.0 GeV 


Hb = 7.5 GeV 


l/R[Ge\] 


^(O)cff 
'-77 


^(o)cfr 


^(O)eff 


^(o)cfr 

'-'8G 


^(O)efr 
^77 


^{0)cff 


200 


-0.236 


-0.076 


-0.192 


-0.053 


-0.169 


-0.040 


250 


-0.264 


-0.100 


-0.223 


-0.079 


-0.201 


-0.068 


300 


-0.282 


-0.114 


-0.242 


-0.096 


-0.221 


-0.086 


400 


-0.301 


-0.131 


-0.264 


-0.114 


-0.244 


-0.105 


SM 


-0.331 


-0.156 


-0.296 


-0.142 


-0.278 


-0.135 



Table 3: Wilson coefficients Cy^ and Cgg in LO for rrit = 167 GeV as functions of 1/R and 
various values of 

We observe a sizable impact of the KK modes on the coefficients Cy^^^^(/X(,) and C^^^{iJ,b) 
although this impact is substantially smaller than on Cj!J{fXw) and C^{fj,w) in (4.9). This 
can be understood by investigating the size of the different terms in (4.6) and (4.7). Setting 
mt = 167 GeV, 1/R = 300 GeV, fib = 5GeV and ai^\Mz) = 0.118 we find 

Cjf *'(//6) = 0.695 c}°Viy) + 0.085 C^^^f^w) - 0.158 (^ly) 

= 0.695 (-0.121) + 0.085 (-0.033) - 0.158 = -0.245 (4.11) 
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to be compared with —0.300 in the SM. In the absence of QCD corrections we would have 
^Tj^^il^b) = c!^^{ijLw) (in that case one has rj = 1) and the KK effects would be very large as 
found in section 3. However, the additive QCD correction present in the last term in (4.11) that 
causes the large QCD enhancement of the B Xsj [36, 37] screens considerably the effects of 
the KK modes. Similar comments apply to C^^^ for which we find 

ciTil^b) = 0.727 ci^'^M- 0.074 Cfifiw) 

= 0.727 (-0.033) - 0.074 = -0.098 (4.12) 

to be compared with —0.144 in the SM. 

We also observe that the strong /Uft-dependence of both coefficients [38, 35] is large. This 
scale-uncertainty in the leading order can be substantially reduced by including NLO corrections. 

4.3 The Branching Ratio 

In the leading logarithmic approximation one has 

where f{z) with z = rric/mh is the phase space factor in Br{B XceVe). 

The corresponding NLO formulae that include also higher order electroweak effects [39] are 
very complicated and can be found in [28]. As reviewed in [30, 31, 32], many groups contributed 
to obtain these NLO results. In our numerical NLO analysis we benefited enormously from the 
computer programs of the authors of [28, 39]. 

In tabic 4 wc show the results for Br{B Xgj) at LO and NLO for the central values of 
the input parameters in [28] and different values of fib- While in the LO the /if, dependence 
is very sizable, it is very small after including NLO corrections. The strong enhancement of 
the branching ratio by additive QCD corrections calculated first in [40] and confirmed in [29] 
is particularly visible at low 1/R for which the Wilson coefficients of the magnetic operators 
are strongly suppressed by the KK modes and the additive QCD corrections become even more 
important than in the SM. This table demonstrates very clearly the importance of higher order 
QCD calculations in the search for new physics. Without them any comparison between the 
SM, the ACD model and the experimental data would be meaningless. 

The strong suppression of the branching ratio by the KK modes is clearly seen. In fig. 3 we 
compare Br{B — > ^^7) in the ACD model with the experimental data and with the expectations 
of the SM. The shaded region represents the data in (4.1) and the upper (lower) dashed horizontal 
line are the central values in the SM for rric/mb = 0.22 {rric/mb = 0.29). The solid lines represent 
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300 400 500 600 700 800 900 1000 
R-^ [GeV] 

Figure 3: The branching ratio for B Xs^ and > 1.6 GeV as a function of 1/R. See text for the 
meaning of various curves. 

the corresponding central values in the ACD model. The theoretical errors, not shown in the 
plot, are for all curves roughly ±10%. 

We observe that in view of the sizable experimental error and considerable parametric un- 
certainties in the theoretical prediction, the strong suppression of Br[B ^sT) by the KK 
modes does not yet provide a powerful lower bound on 1 / R and the values 1/R> 250 GeV are 
fully consistent with the experimental result. It should also be emphasized that Br{B Xg'j) 
depends sensitively on the ratio mc/mi that in (4.2) and table 4 has been set to 0.22 in accor- 
dance with the arguments put forward in [28]. As seen in fig. 3, for a value nic/mb = 0.29 that 
has been used in the past, the branching ratio is smaller by roughly 10% and the lower bound 
on 1/R is shifted above 400 GeV if other uncertainties are neglected. In order to reduce the 
dependence on rric/rab a NNLO calculation is required. Once it is completed and the experi- 
mental uncertainties reduced, Br{B — > Xg^) may provide a very powerful bound on l/i? that 
is substantially stronger than the bounds obtained from the electroweak precision data. 

The suppression of Br{B — > Xg^) in the ACD model has already been found in [3]. Our 
result presented above is consistent with the one obtained by these authors but differs in details 
as only the dominant diagrams have been taken into account in the latter paper and the analysis 
was performed in the LO approximation. 
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= 2.5 GeV 


lib = 5.0 GeV 


Hb = 7.5 GeV 


l/i?[GeV] 


Br(LO) 


Br(NLO) 


Br(LO) 


Br(NLO) 


Br(LO) 


Br(NLO) 


200 


1.54 


2.32 


1.02 


2.30 


0.79 


2.28 


250 


1.92 


2.66 


1.37 


2.65 


1.11 


2.63 


300 


2.18 


2.89 


1.61 


2.88 


1.35 


2.86 


400 


2.49 


3.15 


1.90 


3.15 


1.63 


3.13 


SM 


2.99 


3.57 


2.39 


3.58 


2.11 


3.5G 



Table 4: The branching ratio for B — > in LO and NLO in units of 10 as a function of 
1/i? for iTic/mb = 0.22 and various values of Hb- 

5 B ^ Xs gluon 

5.1 Preliminaries 

The decay 6 — > is governed by the operator Qsg and consequently the value of the coefficient 
CsG is crucial here. In the SM the branching ratio for b ^ sg is strongly enhanced by NLO 
QCD corrections [41] with the result 

Br{b sg) = (5.0 ± 1.0) • 10"^, nie/mb = 0.29 (5.1) 

which is larger by a factor of 2.5 than the LO value. The sizable uncertainty shown in (5.1) 
is due solely to the rcnormalization scale dependence. The other important uncertainty is the 
value of nic/mb- For rric/mb = 0.22 we find using the computer program of [41] 

Br{b sg) = (4.1 ± 0.7) • 10"^, rric/mb = 0.22 . (5.2) 

In fig. 4 we show the results for Br{b sg) in the ACD model as a function of 1/i? for 
three values of Hb and rric/mb = 0.22. In obtaining these results we benefited enormously from 
the computer program of the authors of [41]. As anticipated in section 3, there is a very strong 
suppression of the branching ratio in question by the KK contributions and if the strong fib 
and rric/mb dependences could be put under control and the branching ratio could somehow 
be extracted from the data, the decay B ^ Xg gluon could offer a very powerful constraint on 
1/R. Indeed, even for 1/R = 600 GeV a clear distinction between the SM and the ACD model 
predictions can be made for a fixed /x^. 

Unfortunately, it will take some time before the strong fXb and rric/mb dependences can be 
significantly reduced and the extraction of the branching ratio from the experimental data is 
very difficult if not impossible [41, 42]. Yet, it is not excluded that one day this result could 
become relevant. 
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300 400 500 600 700 800 900 1000 
R-^ [GeV] 

Figure 4: The branching ratio for B Xgg for nic/mb = 0.22 and fib = 2.5, 5, 7.5 GeV. Br{B X^g) 
decreases with increasing ^b- The dashed hnes represent the SM prediction. 

6 B ^ X,ii+ii- 
6.1 Preliminaries 

The inclusive B Xs^i^ decay has been the subject of very intensive theoretical and experi- 
mental studies during the last 15 years. On the experimental side only the BELLE collaboration 
reported the observation of this decay with [43] 

Br{B X,ii+ii-) = (7.9 ± 2.1+2-0) . lo^e _ ^g^^^ 

For the decay to be dominated by perturbative contributions one has to remove cc resonances 
by appropriate kinematical cuts in the dilepton mass spectrum. The SM expectation [44] for 
the low dilepton mass window is given by 

Br{B X,^+/i-)sM = (2.75 ± 0.45) • 10"'' , (6.2) 

where the dilepton mass spectrum has been integrated in the limits suggested by BELLE [43] 
and given in (6.13). 

This cannot be directly compared to the experimental result in (6.1) that is supposed to 
include the contributions from the full dilepton mass spectrum. Fortunately future experimental 
analyses are supposed to give the results corresponding to the low dilepton mass window so that 
a direct comparison between experiment and theory will be possible. 

The most recent reviews summarizing the theoretical status can be found in [32, 44]. The 
purpose of this section is the calculation of the impact of the KK modes on the branching ratio, 
the invariant dilepton mass spectrum and the forward-backward asymmetry. 
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In contrast to the B Xg^ decay the inclusion of NLO corrections to Br{B XsIjl^h~) 
does not require the calculation of QCD corrections to the relevant penguin and box diagrams. 
Detailed explanation of this feature is given in [33, 32, 45]. Consequently in this decay we are in 
the position to present the complete NLO analysis in the ACD model. On the other hand the 
NNLO corrections to B ^ XsH'^n" have been calculated by now within the SM [24, 46, 47] and 
it will be of interest to include these corrections also in the ACD model following the philosophy 
as in the case of NLO corrections to B ^ Xgj and B ^ Xg gluon decays. However, in order to 
see the impact of the KK modes on B ^ XsH'^n~ in a transparent manner we will first present 
a NLO analysis of this decay. 

6.2 Effective Hamiltonian 

The effective Hamiltonian at scales /x = 0{mb) is given by 

Hesib ^ sfi+fi-) = Hes{b ^ S7) - ^^Vtb [Cgv{f^)Q9V + Cioa{Mw)Qioa] , (6.3) 

where we have again neglected the term proportional to V*gVub and 7ies{b 57) is given in 
(4.3). In addition to the operators relevant for B — > Xgj, there are two new operators: 

Qgv = {sb)v-A{P-IJ')v , QlOA = {sb)v-A{P-IJ')A , (6-4) 

where V and A refer to 7^ and 7^75, respectively. They arc generated through the clcctrowcak 
penguin diagrams of fig. 1 and the related box diagrams needed mainly to keep gauge invariance. 
It is convenient to work with the coefficients Cg and Cio defined by 

CMf^) = T^^9(/x), Cioa(m) = ^Cioifi). (6.5) 

Actually Qiqa does not renormalize under QCD and its coefficient is f^ independent with 

~ Y{xt,l/R) 

CioW = (6.6) 

sin Uy] 

and Y{xt, l/R) defined in (3.36) and calculated in [1]. 

The NLO QCD corrections to Cg{fi) in the SM model have been calculated in [48, 49]. 
Generalizing this result to the ACD model we obtain in the NDR scheme 

^r^(M) = Pr"" + ^^^i^ - 4Z{xt, l/R) + PEE{xt, l/R) (6.7) 

sm 6w 

with 

pNDR ^ 2.60 ± 0.25 (6.8) 
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as in the SM. The uncertainty comes from the ni, dependence, with the highest value corre- 
sponding to = 2.5 GeV and the lowest to luj, = 7.5 GeV. The /Xb dependence is substantially 
reduced through the dependence of the perturbatively calculable matrix elements of the op- 
erators involved. Analytic formula for Pq^'-'^ can be found in [48, 49]. The functions Y, Z and E 
are given in section 3. Pe is O(10~^) and consequently the last term in (6.7) can be neglected. 

In table 5 we show the Cg^^{fij,) as a function oil/R for mt = 167 GeV, different values 
of 111, and as{Mz) = 0.118. We observe that the impact of the KK modes on this Wilson 
coefficient is very small due to the approximate cancellation of these contributions to Y and Z. 
Consequently the impact of the KK modes on this decay proceeds almost entirely through the 
enhancement of the coefficient Cio and the suppression of ci^'^j^^ . 



l/R[GeY] 


fib = 2.5 GeV 


Hb = 5.0 GeV 


= 7.5 GeV 


200 


4.476 


4.226 


3.976 


250 


4.434 


4.184 


3.934 


300 


4.413 


4.163 


3.913 


400 


4.393 


4.143 


3.893 


SM 


4.372 


4.122 


3.872 



Table 5: The values of C^^^(//(,) in NLO for mt = 167 GeV as a function of 1/R and various 
values of jib- 



6.3 The Differential Decay Rate 



We are now ready to present the results for the differential decay rate based on the effective 
Hamiltonian in (6.3) and the spectator model for the matrix elements of Qi. 
Introducing 



s = 



z = 



mb 



(6.9) 



and calculating the one-loop matrix elements of Qi using the spectator model in the NDR scheme 
one finds [48, 49] 

2 



d/dsr{b —>■ sfi'^fi ) a 



r(6 ^ cei') 



47r2 



ts 



Vr 



cb 



f{z)K{z) 



U{s) 



(6.10) 



where 



U{s) = {l + 2s) {\Cf{s)\' + |Cio|2) +4 (l + ^) IC-Jf"!' + 124^^^ ReCf (I) (6.11 
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and Cf{s) is a function of s that depends on the Wilson coefficient C^^^ and includes also 
contributions from four quark operators. Explicit formula can be found in [48, 49]. 

Next, f{z) is the phase-space factor for B Xctv and k{z) = 0.88 [50, 51] is the corre- 
sponding QCD correction. In our numerical calculations 0.104. 
Finally, the dilepton mass spectrum is defined as 

^Br{B ^ Xsfi+fM-)is) = Br{B ^ X,ez^)cxp x f{s) = T{s) . (6.12) 

Before presenting the numerical analysis let us note that among the different terms in (6.11), 
|C'|''^(,s)|^ is essentially as in the SM, [Ciol^ is significantly enhanced, while the last two terms 
in (6.11) are suppressed. However, the last term in (6.11) is negative and consequently its 
suppression results in an enhancement of U{s) in addition to the one due to Ciq. 

The formulae presented above served mainly to illustrate the pattern of various contributions 
that is not modified significantly at the NNLO level [24, 46, 47]. However, the latter corrections 
cannot be neglected and we have included them in our analysis. In obtaining these results we 
benefited enormously from the computer program of the authors of [44, 46] that we generalized 
to the ACD model. 




0.01 0.02 0.03 0.04 0.05 0.2 0.4 0.6 0.8 

S s 



(a) (b) 

Figure 5: j~Br{B Xsl+l^) in the SM (dashed line) and the ACD model for R-^ = 250 GeV. (a) 
Low s region and (b) full dilepton mass spectrum. 

In fig. 5 we show T(s) as a function of s for 1/R = 250 GeV. As anticipated before, the 
KK contributions enhance T{s) with respect to the SM, see fig. 5 (b), except for very small 
s where the contribution of the coefficient C^^^^ dominates and its suppression results in the 
suppression of T(s), see fig. 5 (a). 

At this point the following remarks are in order. The theoretical calculations are cleanest 
for So < 0.25 where the NNLO calculations for the inclusive decays have been completed [46, 47] 
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300 400 500 600 700 800 900 1000 
R-^ [GeV] 

Figure 6: Br{B XsfJ,'^^^^) in the SM (dashed hne) and in the ACD model. The integration hmits 
have been chosen as defined in (6.13). 



and resonant effects due to J/ip, ip' , etc. are expected to be small. For this reason as done in 
[44], we will first present the branching ratio for B Xs^~^fJ.~ by integrating over s only in the 
region 



2m, 



< s < 



M 



J/4> 



0.35 GeVy 



(6.13) 



nib J V "^f> 

that gives in the case of the SM [44] the result in (6.2) where the error comes from the variation 
of fib, and mjnib. 

In fig. 6 we show the branching ratio Br[B Xsfi^fJ,^) as a function of 1/R. In obtaining 
these results we followed closely the procedure of the authors of [44, 46] and generalized their 
computer programs to include the KK contributions. We observe a modest enhancement of 
Br{B Xsfi+n') that for l/R = 300 GeV amounts to roughly 12%. 

Now the experimental number in (6.1) refers to the so-called non-resonant branching ratio 
integrated over the entire dilepton invariant mass spectrum. As we do not have any access to the 
experimental analysis and the way the resonance contributions have been removed, we will follow 
the authors of [44] who integrated the theoretical expressions over the whole dilepton invariant 
mass spectrum and obtained in the SM the result Br{B —>■ Xs/i+/i~)sM = (4.15 ± 0.7) • 10^^ . 
For instance for 1/R = 300 GeV we find 



Br{B ^ Xsfi+fi-)ACD = (4.8 ± 0.. 



10" 



(6.14) 



The enhancement of the relevant branching ratio is now stronger than in the previous case 
as the role of Cio becomes more important for larger s. The result in (6.14) is significantly closer 
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to Belle's result in (6.1) than the SM prediction. However, the very large experimental error and 
still sizable theoretical uncertainties in the branching ratio corresponding to the full dilepton 
mass spectrum, also of non-perturbative origin, preclude definite conclusions at present. As we 
stated before it is safer to consider the branching ratio for the low dilepton mass window as 
given in in fig. 6. 

6.4 ForwcLrd-Backwcird Asymmetry 

Of particular interest is the Forward-Backward asymmetry in B ^ Xsfi'^n". It becomes non- 
zero only at the NLO level and is given in this approximation by [52] 

SReC'f(5) + 2C(f'^' 



^fb(s) = — / dcosOi — sgn(cos6'i) = -3Cio 

T{h^ceu)J_i dsdcosUi U[s) 

(6.15) 

with U{s) given in (6.11) and where 0i is the angle between the and B meson momenta in 
the center of mass frame. Similar to the case of exclusive decays [53], the asymmetry A-p^is) 
vanishes at s = sq that in the case of the inclusive decay considered is determined through 

soRe Cf (so) + 20^°^"*^ = 0. (6.16) 

The fact that Afb(s) and the value of sq, being sensitive to short distance physics, are in 
addition subject to only very small non-perturbative uncertainties makes them particularly 
useful quantities to test for physics beyond the SM. 

The formulae (6.15) and (6.16) have recently been generalized to include NNLO corrections 
[46, 47] that turn out to be significant. In particular they shift the NLO value of Sq from 
0.14±0.02 to 0.162ib0.008. In fig. 7 (a) we show the normalized Forward-Backward asymmetry^ 
ApB (s) for the central values of the input parameters that we obtained by means of the formulae 
and the computer program of [44, 46] modified by the KK contributions calculated here. The 
dependence of sq on 1/R is shown in fig. 7 (b) . 

We observe that the value of sq is considerably reduced relative to the SM result obtained 
including NNLO corrections [44, 46, 47]. This decrease originates in the suppression of the 
coefficient c!^^^^ as clearly seen in (6.16). We recall that Cf^ is only weakly affected by the 
KK contributions. For 1/R = 300 GeV we find a value for sq that is very close to the NLO 
prediction of the SM. This result demonstrates very clearly the importance of the higher order 
QCD corrections, in particular in quantities like sq that are theoretically clean. We expect that 
the results in figs. 7 (a) and (b) will play an important role in the tests of the ACD model in 
the fulurc. 



^The normalized FB asymmetry is given as j4fb(s) = r(6 —^ cev) x ^fb(s)/ /^j^ d cos 6i ^^^^^^g'^ ^ 
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(a) (b) 

Figure 7: (a) Normalized Forward-Backward asymmetry in the SM (dashed hne) and ACD for = 250 
GeV. (b) Zero of the forward backward asymmetry ApB the SM (dashed Hne) and the ACD model. 



6.5 Correlation between sq and Br{B — > Xs7) 

In MFV models there exist a number of correlations between different measurable quantities 
that do not depend on specific parameters of a given model. These correlations originate in 
the general property of these models: the new physics enters only through the Inami-Lim (IL) 
functions. As several processes depend on the same IL functions, eliminating these functions 
in favour of measurable quantities allows to derive interesting universal relations between these 
quantities that do not involve parameters specific to a given model. Examples can be found in 
[17, 54]. 

Here we would like to point out a correlation between so and Br{B Xg^) that is present 
in the ACD model and in any MFV model in which new physics contributions to Cg^ and 
its s dependence can be neglected to first approximation. This is the case of the ACD model 
and of a large class of supersymmetric models discussed for instance in [44]. As new physics 
contributions to the Z'^-penguin represented by the function C are suppressed in Cg^ by the 
factor (4 sin^ 0^ — 1), the necessary requirement for the weak dependence of Cg^ on new physics 
is the smallness of these contributions to the B and D functions. This in fact is the case for 
many MFV models. 

If Cg^ is unaffected by new physics contributions and its s-dependence is weak then, as seen 
in (6.16), So is simply proportional to dj'i^^^ with the proportionality factor independent of new 
physics. Consequently sq is proportional to Br{B Xs^), 
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{Br{B X,7) X 10'')2 

Figure 8: Correlation between Br{B Xgj) and sg- Tire straigiit line is a least square fit to a linear 
function. The dots are the resuhs in the ACD model for 1/R = 200, 250, 300, 350, 400, 600 and 1000 GeV 
and the star denotes the SM value. 

and can simply be predicted from the measured value of Br{B — > Xs^) and C^^ given in the 
SM. This correlation is shown in fig. 8. It depends of course on the relation between C^^^°^ and 
Br[B — > Xs7), in particular on mc/nni,. In fig. 8 we have used mc/mi^ = 0.29. The plot in fig. 8 
has been obtained within the ACD model by expressing 1/R in fig. 7 (b) through Br[B — > Xg^) 
by means of the results in fig. 3. The fact that the result of this exercise is a straight line 
confirms our result that the new physics contributions to C^^ and its s dependence are small. 




7 Kl^ TT^e+e- 

7.1 Preliminaries 

There are three contributions to this decay: CP-conserving, indirectly CP-violating and directly 
CP-violating. Unfortunately out of these three contributions only the directly CP-violating one 
can be calculated reliably. In this contribution there are practically no theoretical uncertainties 
related to hadronic matrix elements because {■K^\{sd)v-A\Ki) can be extracted using isospin 
symmetry from the well measured decay — > vr'^e^i^. In what follows we will only consider 
this contribution. Most recent reviews of Ki — > vr'^e'^e" with references to earlier literature can 
be found in [55]. 

7.2 The Branching Ratio 

Generalizing the NLO formula in [56] to the ACD model we obtain 

Br{Ki^ ^ 7rOe+e-)dir = KeilvuXt? {y'^ a + Viv) > (7-18) 
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where ImA* = Im(Ftd^ts), 
VIt{K+) V27r 



/Co — 



6.34 • 10"^ 



and 



yiv = Po + ^^""^^y^^ - V^) + ^^;^(^t, 1/^) , 



sm 



sm 



-Y{xt, l/R) 



(7.19) 

(7.20) 
(7.21) 



with Y{xt, l/R), Z{xt, l/R) and E{xt, l/R) given in section 3. The formula (7.20) has a struc- 
ture identical to (6.7) relevant for B — Xsfi^ ji^ with different numerical values for Pq and Pe 
due to different scales jj, involved. Pe is C(10~^) and consequently the last term in (7.20) can be 
neglected. The next-to-leading QCD corrections to the coefficients above enter only Pq. They 
have been calculated in [56]. One finds Pq = 3.05 it 0.06, where the error comes from the and 
scale uncertainties. Similarly to B — > XsiJ,'^n~, the effect of KK contributions is mainly felt in 

yrA- 

The present experimental bound from KTeV [57] 



Br{Ki_ 



< 3.5 ■ 10"^" {MC.L 



(7.22) 



is still by two orders of magnitude away from the theoretical expectations in the SM and the 
ACD model. 



7.3 Numerical Analysis 

In table 6 we show Br{K\^ 7r*^e^e~)dir in the ACD model as a function of \/R for Pq = 3.05 
and rrit = 167 GeV. We also show there ImA^. To this end we have used the analysis of the 
unitarity triangle of [1]. The enhancement of Br{Ki^ 7r*'e+e~)dir is at most 10% because the 
enhancement of the function Y in y-jA is compensated to a large extent by the suppression of 
ImA(. 



1/P [GeV] 


200 


250 


300 


400 


SM 


ImAt X lO'' 


1.202 


1.250 


1.276 


1.302 


1.333 


Br{Ki^ 7r°e+e-)dir x 10^^ 


4.87 


4.77 


4.69 


4.58 


4.39 



Table 6: Values for ImAt and Br{Ki, 



7rOe+e- 



)dir for different l/R. 
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8 The Ratio e' j e 



8.1 Preliminaries 

The ratio e' je that parametrizes the size of the direct CP violation with respect to the indirect 
CP violation in Ki^ — ^ tttt decays has been the subject of very intensive experimental and 
theoretical studies. On the experimental side the world average based on the recent results from 
NA48 [58] and KTeV [59] and previous results from NA31 and E731 collaborations reads 

e'/e = (16.6 ± 1.6) • lO""^ . (8.1) 

On the other hand, the theoretical estimates of this ratio are subject to very large hadronic 
uncertainties. While several analyses within the SM find results that are compatible with (8.1), 
it is fair to say that the chapter on the theoretical calculations of e' je is certainly open. Most 
recent reviews with references to original papers can be found in [60, 61, 62]. 

In view of this situation, our strategy will be to choose various sets of the values of the main 
non-perturbative parameters B^^"^^ and B^^'^^ for a fixed strange quark mass ms(mc) and to 
investigate for each chosen set (-Bg^^^^ , B^^'^^ ) the impact of the KK modes on e'/ e. 

8.2 Basic Formula 

The formula for e'/e of [63, 64] generalizes to the ACD model as follows: 

^ = lmXfF,,{xt,l/R) , (8.2) 

where 

F,,ixt, l/R) =Po + Px X{xt, l/R) + Py Y{xu l/R) + Pz Z{xt, l/R) + Pe E{xu l/R) (8.3) 

with the functions X, Y, Z, E given in section 3. Strictly speaking the functions X and Y 
entering (8.3) are given hy X = C + and y = C + B^'^'^) and consequently they differ 

slightly from the ones defined in (3.35) and (3.36) that involve B'^^ and i?^^ instead of 
and B^'^'^\ respectively. However, this difference caused only by the KK contributions to the 
box diagrams is fully negligible and the numerical values for X and Y given in table 1 can also 
be used here. 

The coefficients Pi arc given in terms of the non-perturbative parameters B^^"^^ and B^^"^^ 
and the strange quark mass ms(mc) as follows: 

P, = rf + rf i?6 + rf)i?8 , (8.4) 
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where 



^6 = 



(1/2) 



121 MeV 



rrisimc) +md(mc)_ 



i?8 = B. 



(3/2) 



121 MeV 



msirric) +md(mc) 



(8.5) 



(1/2) (3/2) 

Bq and parameterize the matrix elements of the dominant QCD-penguin (Qg) ^-iid the 

dominant electroweak penguin (Qs) operator, respectively. The numerical values of r^^\ rf'^ 
and rf'^ for different values of at /i 
in the NDR renormalization scheme in table 7. 



nic have recently been updated in [62] and are given 





M 


4 = 310 MeV 


M 


4 = 340 MeV 


M 


4 = 370 MeV 


i 


i 


i 


i 


i 


i 


i 


i 


i 


^(8) 
i 



X 
Y 
Z 
E 


-3.167 
0.546 
0.395 
0.626 
0.189 


14.781 
0.026 
0.107 
-0.020 
-1.705 


1.815 



-12.574 
0.552 


-3.192 
0.537 
0.385 
0.673 
0.179 


15.974 
0.029 
0.113 
-0.021 
-1.801 


1.686 



-13.226 
0.592 


-3.215 
0.529 
0.376 
0.723 
0.169 


17.277 
0.032 
0.120 
-0.023 
-1.903 


1.546 



-13.927 
0.634 



Table 7: Coefficients in the formula (8.4) for various in the NDR scheme [62]. 



8.3 Numerical Analysis 

The values of ImAt for various values of 1/i? are given in table 6. In table 8 we show the 
corresponding results for e'/e as functions of the non-perturbative parameters in question. To 
this end we have set ms{mc) = 100 MeV which corresponds to m5(2GeV) = 85 MeV, in the 
ballpark of some most recent lattice calculations with dynamical fermions [65]. The last year 
lattice values [66] and the QCD sum rules values [67] are slightly higher: ms{mc) = 115ib20MeV. 

The results in table 8 are self explanatory. Dependently on the values of the parameters 
Sj^/^^ and one finds the value oi 1/R that is favoured by the experimental data. Only 

for the case {Bq'^^\b^^'^^) = (1.0,1.0) the results for e'/e are outside the experimental range 
in (8.1) independently of 1/R considered. As for 310 MeV < A^ < 370 MeV the ratio e'/e is 
approximately proportional to A^, the results in table 8 can be at most changed by ±10%. 

(1/2) (3/2) 

Clearly, the hadronic uncertainties in {Bq ,B^' ') and rus {rric) preclude any determination 
of 1/i? from the data on e'/e at present. As generally e'/e is suppressed by the KK mode 
contributions relatively to the SM, the ACD model is disfavoured by e'/e unless {B^^'^\ B^^'^^) 
and ms{mc) or generally the matrix elements of Qq and Qs are such that the SM estimates are 
above the data. This is the case for sufficiently low values of ms{mc) and sufficiently high values 



28 





o(l/2) 


= 1.00 


o(l/2) 


= 1.15 


r(1/2) 


= 1.30 


1/R 


Bf^^ = 0.8 


Bf/^^ = 1.0 


= 0.8 


5f ^) = 1.0 


= 0.8 


= 1.0 


200 


9.1 


5.6 


12.7 


9.2 


16.4 


12.9 


250 


11.0 


7.7 


14.8 


11.5 


18.6 


15.3 


300 


12.1 


9.1 


16.0 


12.9 


19.9 


16.8 


400 


13.4 


10.6 


17.3 


14.5 


21.3 


18.5 


SM 


15.2 


12.8 


19.. 3 


1G.8 


23.1 


20.9 



Table 8: The ratio e'/s in units of 10"^ for rrit = 167 GeV, = 340 MeV, ms(mc) = 100 MeV 
and md{mc) = 6 MeV as a function of 1/i? in GeV and various values of B^^'^^ and B^^^\ 

ofB^^^^\ For these cases a suppression of e'/e by the KK modes found here is a welcome feature 
of the ACD model. 



9 Summary and Outlook 

In this paper we have calculated for the first time the contributions of the Kaluza-Klein (KK) 
modes to B ^ Xg gluon, B — Xg^'^fi" and — 7r°e+e~ and to the CP- violating ratio e'/e 
in the Appelquist, Cheng and Dobrescu (ACD) model with one universal extra dimension. We 
have also analyzed B — Xgj that has been considered in the past [3]. While our result for 
this decay differs in details from the latter paper, we confirm the suppression of the relevant 
branching ratio found by these authors. As a byproduct we have given a list of the required 
Feynman rules involving photons and gluons that have not been presented in the literature so 
far. Moreover we have generalized the background field method to five dimensions. 

The nice property of this extension of the SM is the presence of only a single new parameter, 
1/R. This economy in new parameters should be contrasted with supersymmetric theories and 
models with an extended Higgs sector. Taking 1/R = 300 GeV our findings are as follows: 

• The short distance function Z, relevant for B — Kl — Tr^e'^e" and e'/e is 
enhanced by 23% relative to the SM value. 

• The functions D' and E', relevant for B ^ Xgj and B Xg gluon are suppressed by 36% 
and 66%, respectively. The corresponding effects in the function D are negligible and in 
E phenomenologically irrelevant. 

• Br{B Xg^) is suppressed by 20%. The phenomenological implications of this result 
depend sensitively on the value of rric/mh and on the experimental data, with the lower 
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bound on 1/R being stronger for mc/m6=0.29 than for mc/mi,=0.22. In fact, in the latter 
case the suppression of Br{B Xg^) could be welcome and an upper bound on l/i? could 
in principle be found when the experimental and theoretical uncertainties decrease. 

• Br{B Xs gluon) is suppressed by 40%. The phenomenological relevance of this result 
is, in view of large hadronic uncertainties and the difficulty in extracting this branching 
ratio from the data, unclear at present. 

• The perturbative part of Br{B is enhanced by 12%. This could be a welcome 
feature of the ACD model in view of the data from Belle, that is above the SM expectations. 
However, more interesting is the shift of the zero in the A^b asymmetry from sq = 0.162 
to So = 0.142, that in view of small theoretical uncertainties could turn out to be a very 
important test of the ACD model. 

• We have pointed out a correlation between the zero sq in the ^fb asymmetry and Br{B 
Xgj) that should be valid in most models with minimal flavour violation. This correlation 
is shown in fig. 8. 

• Br{KL — > 7r°e"'"e~) is enhanced by less than 10%. 

• e' / £ is suppressed relative to the SM expectations with the size of the suppression depend- 
ing sensitively on the hadronic matrix elements. Taking ms(2 GeV) = 85 MeV as suggested 
by the most recent lattice calculations, we find that for several sets of {Bq'^'^\b^^'^^) the 
ACD model can be made compatible with the experimental data, see table 8 for details. 

These findings should be compared with the ones of [1], where for 1/R = 300 GeV the 
following enhancements relative to the SM predictions due to enhanced Z°-penguins have been 
found: K+ tt+uU (9%), i^L ^ 7r°z/i^ (10%), B ^ Xdvv (12%), B Xsuv (21%), Kl 
fijl (20%), Bd ixjl (23%) and Bs nfi (33%). 

We observe that the impact of the KK modes on B ^ Xgj is comparable to the one on 
B XgVV, Kl up. and B^, nJi. The corresponding impact onB ^ Xg gluon is significantly 
larger and comparable to the one on Bg M^A*"- On the whole, the main effects of the KK 
modes are felt in 2^°-penguins, 7-magnetic penguins and chromomagnetic penguins. 

Moreover, the low compactification scale 200 GeV that was fully compatible with the data 
on the decays considered there, is excluded by the B Xgj decay. For 1/R = 200 GeV the 
branching ratio Br{B Xgj) is suppressed relative to the SM by a factor of 1.6. This is clearly 
excluded. Of considerable interest is an even stronger suppression of B ^ Xg gluon. 

Interestingly, the present data on if + tt^uu and B Xgfi'^iJ,~ are somewhat above the 
SM expectations and the enhancement of the .Z^°-penguins by the KK modes could be welcome 
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in this context. On the other hand, the present central experimental values for Br{B — Xgj) 
are somewhat below the SM expectations and the suppression of the magnetic penguins by the 
KK modes could also be welcome. However, in all these cases the experimental and theoretical 
uncertainties have to be reduced in order to be definitive about the need for these enhancements 
and suppressions. 

Whether the suppression of B ^ Xg gluon and of e'/e by the KK modes relative to the SM 
expectations is a welcome feature depends very strongly on hadronic uncertainties and in the 
case of B ^ Xg gluon on the data. 

Possibly, the most interesting result of our paper is the sizable downward shift of the zero 
(so) in the ^fb asymmetry. It should be emphasized that this shift has a definite sign and the 
theoretical uncertainties in its value are small. 

Our present analysis combined with our previous paper [1] completes the study of the most 
interesting FCNC processes in the ACD model. With only one additional parameter, the com- 
pactification scale, the pattern of various enhancements and suppressions relative to the SM 
expectations could be uniquely determined: 

• Enhancements: Kl — ^ ir^e'^e' , AMg, — tt'^vv, Ki, — > tt^uu, B — ^ X^yv-, B — > Xsvu, 
Kl — ^ lJ''^fJ''' , Bfi — jJL^ n~ ., B — ^ XsjJL^ jJi" and Bg — /J-'^/J-" ■ 

• Suppressions: B — > Xsj, B Xg gluon, the value of sq in the forward-backward asym- 
metry and e' /e. 

Whether these enhancements and suppressions are required by the data or whether they exclude 
the ACD model with a low compactification scale, will depend on the precision of the forthcoming 
experiments and the efforts to decrease the theoretical uncertainties. 
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A Background Field method in 5 dimensions 



In the calculation of the off-shell amplitude, we have used the method of background fields. 
As in 1-PI diagrams the background fields only appear on the external legs, there is no need 
to fix the gauge for the background gauge fields. With the appropriate choice of the gauge 
fixing for the quantum gauge fields, we can ensure the invar iancc of the effective action with 
respect to background field (BF) gauge transformations. The general procedure in 4 dimensions 
is described in [25] . Here we only state the SM case and modify the gauge fixing to fit the ACD 
model. 

The starting point is the action S[ip] before gauge fixing, where ip denotes all gauge and 
matter fields in the action. We get the BF action by the transformation 

S[^]-^S[^I; + ^], (A.1) 

where we have introduced the background fields -0. It is modified by the gauge fixing to 

SBF[^,i']=S['iP + i;]-^J d'^x gg + ghost terms, (A.2) 

where gg stands for all gauge fixing functionals. As the terms with ghost fields are not relevant 
for us we will not consider them here. 

In the 4-dimensional electroweak SM, a convenient choice for the gauge fixing functionals is 

ObMB, 0, 0] = d^B^' - ic^i\ (<^t0 _ ^t^) ^ (A.3) 

^1,SM k <!>, ^] = d^A-f^ + 526""'=^^^-'^ - i92^l - </.t(7>) , (A.4) 

with the Higgs fields 

4>=l=(v + i; + ixV) (^^, = 71 + (i) ' (^-^^ 

where cr" are the common Pauli matrices. According to (A.2), the gauge fixing part of the 
Lagrangian is 

jo.GF = -^gBgB-^gm. (a.g) 

With this specific gauge fixing, the BF action is invariant under the BF gauge transformation^ 
<^BF^^ = ^d^P, (A.7) 
fepi^ = r'^Ay + -d^a^ (A.8) 

92 

6bf4> = i (a" Y + ^P) 4>, (A.9) 



^We omit the BF gauge transformation of the fermions which is just an ordinary gauge transformation. 
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combined with a transformation of the quantum fields 



M» = r^'^Ay, (A.io) 

,5</. = z(a'^^ + 1/3^0, (A.11) 

where (A.IO) and (A. 11) are just unitary rotations in the functional integral. 
The analogous BF gauge transformation in 5 dimensions is 

SbfBm = ^dnP, (A.12) 
SbfAI, = r'^Ai.a'' + ^dMa-, (A.13) 

92 

Sbf4> = i f Y + 4> (A.14) 



combined with 



5AIj = f^'Aijo^, (A.15) 
S^ = i(a-^ + U\^. (A.16) 



2 2' ^ 

Compared to 4 dimensions, the couplings g' and §2 have been replaced by their 5-dimensional 
analogs g' and g2- All fields and a" and (3 are now also functions of the extra coordinate y and 
the vector index M can take the values Af = 0, 1, 2, 3, 5. 

In the ACD model, we have to add —^d^B^ to (A. 3) and —^95^5 to (A. 4) in order to 
diagonalizc the mass matrices of the bosonic modes [1]. However, this spoils BF gauge invariance, 
so we add another term to (A. 4) to fix this and get 

Gb,acb[B, (t>, 4>] = QbMB, ^] - ^d^B^, (A.17) 

^.4,acd[^, ^, 0, 4>] = GaMA a, 0, 4>] - ^d^At - C92e''''AlAl (A.18) 

where it is understood that g' and g2 have been replaced by g' and 52- 

The use of (A.17) and (A.18) ensures invariance of all 1-PI diagrams under BF gauge trans- 
formations in 5 dimensions. However, as we are only interested in external zero-mode fields, the 
last term in (A.18) does not affect our calculation. 

In the diagrams for the calculation of the functions E and E', the external gluon couples 
only to the quarks. As there are no fermions involved in the gauge fixing, the background gluon 
couples to quarks just like a quantum gluon, and there is no need to specify a gauge fixing 
functional for QCD. For completeness, we note that 

Gg,acd[G, G] = d^G'^'^ + gsr^'C'^C^^ - ^d^G^ - ^gsf'^GlGi (A.19) 
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would be an appropriate choice. Here / are the SU{3)c structure constants and Qs is the 
strong coupUng constant. 



B Feynman Rules in the ACD Model: Photon and Gluons 

In this section we hst the Feynman rules needed for the calculations in this paper except for 



described in appendix A. The rules for vertices involving only quantum fields are the same as 
in the conventional 5d i?^-gauge described in [1]. 

In order to simplify the notation, we omit the KK indices of the fields. There is no ambiguity 
because in one-loop calculations at least one field is always a zero-mode. In the vertex rules 
given below, this is the photon A, the gluon G and their background equivalents. Due to KK 
parity conservation, the other two fields have equal KK mode number, i.e. either zero or n > 1. 

Fermion zero-modes have substantially different Feynman rules than their KK excitations. 
The up-type quarks Qi and Ui are always supposed to be (n > l)-modes, while the zero-mode 
is labeled u. The generation index i can take the values i = u,c, t. 

In the vertices below, stands for the scalar modes and a^. All momenta and fields 
are assumed to be incoming. The Feynman rules for the vertices are: 



those already given in [1] . The Feynman rules are derived in the 5d background field ii!^-gauge 



5^ 




AW^G 



C=l 



AW-G+ : C = -l, 



(B.l) 



AW+a" : C = 0, 



AW'a 



.+ 



C = 0, 



(B.2) 



AW+G- : C = 0, 
Aw^a- : C = 0, 



AW-G+ : C = 0, 
Aw-a+ : C = 0. 



(B.3) 
(B.4) 
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= -ig2Sw{k2 - ki)nC. 



AG+G- : C = l, 
AG+a- : C = 0, 



Aa+a- : C = 1, 
Aa+G- : C = 0, 



and the same values of C for the analogous vertices with a background photon A. 



























AW+W- : 


C = 


9nu{k2 -ki)x 


AW+W- : 


c = 


9ij,u{k2 -ki + 




+ 9\u{k3 - k2)^. 


F2 












= ig2SwiiJ,c. 











AuiUi 

AQiQi 

AQiUi 



2 

3' 
2 

3' 

C = 0, 



c 
c 



AUiUi : C 



3' 



(B.5) 
(B.6) 



(B.7) 
(B.8) 
(B.9) 



AUiQi-. C = 0, 



(B.IO) 
(B.ll) 
(B.12) 



and the same values of C for the analogous vertices with a background photon A. 
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GuiUi : C = l, (B.13) 

GQiQi : C = l, GUiUi : C = 1, (B.14) 

GQiUi : C = 0, GUiQi : C = 0, (B.15) 

and the same values of C for the analogous vertices with a background gluon G. 
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